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Abstract 


This paper presents a decentralized control design method for large complex flexible 
structures by using the idea of joint decoupling. The derivation is based on a coupled 
substructure state-space model, which is obtained from enforcing conditions of interface 
compatibility and equilibrium to the substructure state-space models. It is shown that by 
restricting the control law to be localized state feedback and by setting the joint actuator 
input commands to decouple joint degrees-of-freedom from interior degrees-of-freedom 
the global structure control design problem can be decomposed into several substructure 
control design problems. The substructure control gains and substructure observers are 
designed based on modified substructure state-space models. The controllers produced 
by the proposed method can operate successfully at the individual substructure level as 
well as at the global structure level. Therefore, not only control design but also control 
implementation is decentralized. Stability and performance requirement of the closed- 
loop system can be achieved by using any existing state feedback control design method. 

A two-component mass-spring-damper system and a three-truss structure are used as 
examples to demonstrate the proposed method. 


I* Introduction 

Control of flexible structures has gained much research interest since the space shut- 
tle transportation system became reality and construction of large structures in space is 
no longer a dream. Many structural control methods and algorithms have been proposed 
in the past decade. However, the ap plication of decentralized control to flexible structures 
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has not been pursued extensively. A large space structure must be built incrementally 
with components shipped into space sequentially during several shuttle missions. Also, 
for operational purposes, space structures may need to be connected and disconnected 
routinely in space. Therefore, it is desirable to have controllers that can function well 
on each individual structure and also on the interconnected structures. Decentralized 
controllers designed based on structure components will fit to meet such purpose. 

Most of existing decentralized control designed methods were developed for ap- 
plication to electrical engineering systems or economics systems, rather than to flexible 
structural dynamics systems (see, for example, Refs. [1-3]). In general, structural dy- 
namics systems are strongly coupled (in physical coordinates), which makes it difficult to 
design and to implement decentralized control. Existing decentralized control methods 
for flexible structures either adopt or extend the concepts and methodologies developed 
for other type of systems. For instance, the method in Ref. [4] extends the concept of 
decentralized fixed modes of Ref. [1] to structural control problem. The decentralization 
of the control problem is based on modes instead of physical components of the structure. 
The controlled component synthesis (CCS) method of Ref. [5] adopts the concept of over- 
lapping decomposition in Ref. [3] to design controllers for physical structural components. 
Because overlapping decomposition requires information about adjacent components, the 
CCS method is not a strict decentralized approach. The substructure controller synthe- 
sis (SCS) method in Ref. [6] is a decentralized control design method which uses natural 
decomposition of structural dynamics systems. Both CCS and SCS methods employ 
interface compatibility conditions to assemble the substructure controllers into a global 
controller for the assembled structure. However, because interface equilibrium conditions 
are not considered in the assembly process, the global controller does not guarantee sta- 
bility of the closed-loop system, which is a major disadvantage of both CCS and SCS 
methods. 

This paper presents a decentralized control design method based on the idea of 
joint decoupling. The method is developed specifically for decentralized control of large 
complex flexible structures. The structure to be controlled is first decomposed into sev- 
eral substructures, which can be physical components that are to be connected together 
to form the assembled structure, or they can simply be the result of applying imagi- 
nary” cuts to the structure. It is assumed that there are collocated actuators and sensors 
at every interface (or joint) degree- of- freedom. Then, by enforcing conditions of com- 
patibility and equilibrium to the joint inputs and joint outputs, state-space models of 
substructures are combined together to form a coupled substructure state space model. 
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Then, based on the coupled substructure state-space model it is shown that control de- 
sign can be decentralized by restricting the control gain matrix and the observer gain 
matrix to have a specific form. The whole design concept can be described as follows. 
The control commands for the interior actuators (i.e., actuators that are not located at 
the joints) of each substructure are chosen to be localized state feedback. Then, the 
joint actuator commands are set to cancel out all the forces acting on the joint degrees- 
of-freedom which are induced by interior actuator inputs and by vibration of interior 
degrees-of-freedom. The result is that the joint degrees-of-freedom are decoupled from 
the interior degrees-of-freedom and, therefore, there is no interaction force transmitted 
through the joints from one substructure to another substructure when the substructures 
are connected together. The controller designed by using the idea of joint decoupling is 
called the Joint Decoupling Controller. 

The proposed Joint Decoupling Controller has several advantages. First, the de- 
sign procedure is completely decentralized. It requires no information about the other 
substructures in order to do control design for one substructure. Second, the control 
implementation is also decentralized. The controllers can operate successfully at the in- 
dividual substructure level before assembly and also at the global structure level after 
assembly. This feature makes the proposed method extremely useful for active control of 
space structures that need to be connected and disconnected on a routine basis for op- 
erational purposes. No controller redesign or controller shut-off-and-turn-on is necessary 
before, during, or after the connection process. Third, the proposed method divides a 
large-scale control problem into several small-scale subproblems, and, hence, computa- 
tionally it is more efficient than a centralized control design approach. The control law 
can be determined by using any state feedback control design method, e.g., the LQR 
method or pole placement method. 

This paper is organized as follows. The coupled substructure state-space model is 
derived in Section II. The design procedure for a decentralized control gain matrix and the 
design procedure for a decentralized observer system are presented in Sections III and IV, 
respectively. Also included in Section IV is a summary of design steps for the proposed 
Joint Decoupling Controller. A six degrees-of-freedom mass-spring-damper system and 

a three-truss structure are used in Section V as examples to demonstrate the proposed 
method. 
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II. The Coupled Substructure State-Space Model 


In this section, we will derive the coupled substructure state-space model for flexible 
structures that are composed of many substructures. Without loss of generality, we will 
consider the two-component structure shown in Fig. 1. Figure 1(a) shows the assembled 



Figure 1: A two-component structure. 

structure as an entity, which also will be referred to as the global structure. Figure 1(b) 
shows the decoupled two substructures a and /?. These two substructures can be physical 
components that are to be connected together to form the structure, or they can be the 
result of applying an “imaginary cut” to the structure. Locations of actuators and sensors 
are denoted by u and y, which are divided into two groups. Those actuators and sensors 
located at the joint (or interface, where the substructures are to be connected together) 
are denoted with a superscript J. Those actuators and sensors located at unconnected 
points (or located at “interior” points of each substructure) are denoted with a superscript 
I. In order to describe interface compatibility and equilibrium conditions in terms of 
input and output vectors, it is assumed that there are a pair of collocated actuator and 
sensor at every joint degree- of-freedom. Also, the output measurements at the joint 
coordinates are assumed to be accelerations for a special purpose to be discussed later. 
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Let the state-space models of a-substructure and /^-substructure be described, 
respectively, by 




x a + 


D'J D { 3 

Di' D J J 


if = Apxg + [ Bj, ] | ”j j 
f*i LfCjl JD' 0 ' D'S 1 

lv«l 1 ci J . d " 


These two state-space models can be derived from finite-element modeling of the two 
substructures. For instance, let the finite-element model of a-substructure be described 


(m 


( vl \ = [ H l 0 ] r wi 1 

1 vi I o / J \ wi } 


where w a and are the interior and the joint displacement coordinates, respectively. 
M a , Z a , and I< a are the substructure mass, damping, and stiffness matrices. P £ is 
the interior actuator distribution matrices. H ! Q is the interior accelerometer distribution 
matrices. (For more general cases, the interior outputs may also include displacements 
and velocities.) The identity matrix I in the input and output matrices indicates that 
there are actuators and sensors collocated at every joint degrees-of-freedom. In first-order 
state-space form, the system in Eq. (3) can be represented by 


-M~ x I< a -M~'Z a 


iih 


f 51 1-*.-*. f ^ ® 


where the state vector is defined by 
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For convenience of derivation in this paper, the substructure state vector is rearranged 


as 
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and the system matrices are rearranged accordingly. Other than finite-element model- 
ing, the substructure state-space models also can be realized directly from substructure 
experimental data by using system identification tools. 

The two state-space models in Eqs. (1) and (2) describe the dynamics of the two 
substructures when they are completely decoupled. When the two substructures are 
joined together to form the global structure, the input and output vectors at the interface 
of the two substructures are no longer independent. They must satisfy two conditions: 
compatibility and equilibrium. The output vectors at the interface are constrained by 
the compatibility equation: 

( 6 ) 


y J a 


J J 

yp = y 


which says that the physical motion of the two substructure at the interface must be the 
same. The input vectors at the interface are related by the equilibrium equation: 


U J = U J a + Up 


(7) 


which says that the sum of the internal forces at the interface must be equal to the 
external force applied there. The two substructure state-space models can be coupled by 
enforcing interface compatibility and equilibrium conditions. First, rewrite the bottom 
part of the output equations in Eqs. (1) and (2) as 


u J s = {Dj 1 )"' (yj - Cjxg - 

Adding the above two equations and applying the equilibrium and compatibility condi- 
tions in Eq. (7), we obtain 


u J = [K J )-‘+ (DiS-r 1 ] V J - {Dt j yX C i*° + Di'O - {Di J )-\cix 0 + D}' «Q 

The above equation can be rewritten as 

y J =Dt J S-'D J / k + « J f,+ (Di j y'(Cix a +Di'u' a ) + (D J /) (9) 
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where S = (. D J a J -f D 0 J ^j and u 3 is replaced by u J a + u J 0 . Finally, by substituting Eq. (9) 
into Eqs. (8) and then substituting the results into the two substructure state-space 
models in Eqs. (1) and (2), we obtain the following “combined” state-space model for 
the two substructures: 

A a - B J a S~'C J a B J a S-'Ci 


{i;l- 


B J „S-'Cl 


A„ - BjS 


'■ 'i 1 i 
~ 1( ^0 . I xp / 
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b}s-'d " 
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b j ,s-'d j „ j \b'-b1s-'d j ' 


a 1 


( ' 


■Ci-D^ J S-'C J a D”S-'C} ■ 

yt 


D J „ J S-'Ci Dl J S~'C l 

u 


Dj J S-'Ct Di J S-'C J e 

l y r o , 


. D'/S-'Ci Cl-D'JS-'C 1 ,,. 



{:;} 


(10) 
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1 

< 
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i 

4 

. D'/S-'Di 1 
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In short notation, Eq. (10) will be represented by 

x = Ax + Bu 
y = Cx + Du 


(11) 


where A, B, C, D, x, u, and y are clearly defined. 

The state-space model in Eq. (10) will be referred to as the coupled substructure 
state-space model. It describes the dynamics of the two substructures in Fig. 1(b) with 
compatibility and equilibrium conditions satisfied at the interface. Therefore, it can be 
considered as a state-space representation of the global structure in Fig. 1(a). In fact, 
there is a relationship between the coupled substructure state-space model and the state- 
space model of the global structure. Let the state-space model of the global structure in 
Fig. 1(a) be described by 

x — Ax -f Bu 

y = Cx + Du ^ 
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with the state, input, and output vectors defined by 



x 1 1 

f U a ' 
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We see that x and x are related by a linear transformation 
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x = Tx 


( 13 ) 


(14) 


where T will be referred to as the state coupling matrix. 


Also, u and u are related by 
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y and y are related by 
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O 
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1 
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0/0 
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I 

o 

0 
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T\ and T 2 will be referred to as the input coupling matrix and the output coupling matrix , 
respectively. It should be noted that in general T\ / 7^, although it appears to be, 
because the number of interior sensors and the number of interior actuators can be 
different. By using the state, input, and output coupling matrices, it can be shown that 
the global structure state-space model and the coupled substructure state-space model 
are related by 


TA = AT , TBTi = B , T 2 C = CT , T 2 DT, = D (17) 


In the preceding derivation, we have assumed that both D J a 3 and Dp J are invertible. 
This is true when the outputs at the interface are accelerations. As can be seen from 
Eqs. (3) and (4), when the interface output is y 3 — w 3 we get 


D" = [ 0 I ] AC 1 


0 
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which is invertible because M a is positive definite. (For most structures, the mass ma- 
trix is positive definite.) Since accelerometers are the most frequently used sensors in 
practical situations, this assumption is not a serious restriction. Even when the sensing 
devices used at the interface are displacement sensors or velocity sensors instead of ac- 
celerometers, it is still possible to perform coupling of structure state-space models by 
making modifications to the substructure state-space models. 

In the following sections, a method for designing decentralized controllers for the 
global structure will be derived. It will be shown that the controller gains for joint 
inputs and the observer gains for joint outputs can be chosen such that the coupling 
between substructures is eliminated. By doing so, the global structure control problem 
is decomposed into two substructure control problems. 

III. The Design of Decentralized Controller Gain Matrix 


According to the separation principle, controller gain matrix and observer gain 
matrix can be designed separately. In this section, we will derive the procedure for 
designing decentralized (or localized) state feedback gain matrix for the global structure. 
As discussed previously, the coupled substructure state-space model in Eq. (10) is a 
state-space representation of the global structure. Therefore, the whole derivation will 
be based on the coupled substructure state-space model. 

The key idea behind the decentralization method presented here is to use the joint 
actuators to eliminate the interaction between substructures. Let the joint actuator 
commands be chosen as 


< = - {Di J Y' {Ci + Di’ u' a ) 

u i — ~ (^j J ) {ps + Bp u p) 


After algebraic manipulations, substitution of the above into the state equation in Eq. (10) 
yields 


J 1 _ Aam 0 ( X a 1 Bl m 0 f U* a ) 

1 I 0 A9m J 1 X0 j 0 B r 0m J \ U { 0 j 


( 19 ) 


where 


= A„ -Bi (DP)'' Cl 

A pm — Ap — Bp {p J p J ^ Cp 


BL = fl' - Bi (dp)-' DP 
Bl m = B‘ - Bl (DP)'' D" 
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Subscript m is used to denote that the above matrices are modified substructure matrices. 
Note that the system in Eq. (19) is decoupled. This decoupling of the closed-loop state 
equation suggests a decentralized control design procedure, which can be described as 
follows. First, restrict the feedback commands for the interior actuators to be localized 
state feedback as described by 

U a ~ G l a X Q , U0 = GpX0 (21) 


The gain matrices G ! a and Gjj are to be designed based on the following modified sub- 
structure state- space models : 


%a — A am ^a "I" B am U a 
x 0 = A 0m x 0 + Bp m ujj 


( 22 ) 


Then, according to Eq. (18), set the feedback commands for joint actuators to be 


<=- (Cl + Dl'G'J) *. = Gl x. 

<4 = - ( D e J ) 1 ( c i + D i' G h) x t> = G i x t> 


Finally, substitute this set of actuator commands in Eqs. (21) and (23) into Eq. (15) 
and use the input coupling matrix to obtain a global feedback gain matrix for control 
implementation on the global structure. 

For the decentralized control design approach described above to be successful, 
both modified substructure systems (A am , B r am ) and (Ap m , Bp m ) must be controllable. 
Or, at least, all the unstable poles of A arn and Ap m need to be controllable so that the 
modified substructure systems are stabilizable. However, this is not true. Substitution of 
the joint actuator commands in Eq. (18) into Eq. (8) gives = 0 and yp = Wp = 0, 

which means that all the joint degrees-of- freedom are free of forces. Physically, the joint 
actuator commands in Eq. (18) cancel out all the internal forces acting on the joint 
degrees-of- freedom. The internal forces acting on the joint degrees-of-freedom are induced 
by interior inputs and vibration of the substructures. Cancellation of internal forces 
decouples the dynamics of the joint degrees-of-freedom from that of the interior degrees- 
of-freedom. By doing so, there will be no interaction between substructures when the 
two substructures are connected. The only problem is that the joint degrees-of-freedom, 
even though free of forces, still can have rigid body motion. This rigid-body motions of 
the joint-degrees-of-freedom are not controllable by interior actuators, since all control 
forces coming from interior actuators to the joint degrees-of-freedom will be cancelled 
by the joint actuator commands in Eq. (18). As a result, the substructure closed-loop 
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systems and the global structure closed-loop system attain a pair of unstable double poles 
at the origin for every joint degree-of-freedom. A simple approach to control the rigid 
body motions of joint degrees-of-freedom is to introduce an augmented joint actuator 
command. Details about the design of interior actuator feedback gain matrices, G J a and 
Gp, and the augmented joint actuator command are discussed in the following. 

In physical coordinates, the modified state equation for ar-substructure can be 
partitioned as 



A 11 A IJ 1 ( x 1 1 B n 

am ™am I x a I D am / 

A JI A JJ | J f + dJI u c y 

am ^orm J t x a ) G am J 


(24) 


in which the state vector x a is partitioned according to Eq. (5). As discussed previously, 
the uncontrollable state in Eq. (24) is x J a . Therefore, (A" m , Bj/J is a completely control- 
lable system, presuming that the given interior actuator setting makes the substructure 
controllable. We can set 



with G r J to be designed based on the subsystem x * = A^x^ + B^u^. Any existing 
state feedback design method, e.g., the LQR control theory or pole placement method, 
can be used to determine G r J. Similarly, we can set Gf, = [ Gff 0 ] and let GjJ be 
designed based on ip = Ap^ip + Bj^Up. By doing so, the controllable poles of A am and 
Ap m are stabilized or controlled by u ' and u J. If the substructure state-space models are 
identified from experimental data rather than derived from finite-element models, the 
system equations are in general not in physical coordinates. Then, there is no way to 
tell which part of the system matrices are associated with the interior degrees-of-freedom 
or with the joint degrees-of-freedom. In this case, it is necessary to use existing system 
realization algorithms to separate the system equations into controllable subsystems and 
uncontrollable subsystems. Then, G l J and Gjf can be designed based on the controllable 
subsystems. 

To control the rigid-body motions of the joint degrees-of-freedom, an augmented 
joint actuator command is introduced. Let 


u t = Gix a + ui 

Up — GpXp + u J p ^ ^ 

where u a and Up are the augmented joint inputs. Substitution of the above into Eq. (8) 
yields 
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*i = {d j /y' vi = { d j /y' u> j 0 


We can set the augmented joint inputs to be 

*t = (ot J ) "[-[«?] -[*oa ] ] { % } - *i 

“fl = ( D « J ) [ - [ n i ] - [ 2C,nj 1 ] { wj } - x i 


( 27 ) 


where J 2^,0, ] and [ 0? ] are diagonal matrices with the design parameters £,• and fl,, 
i = 1,2, ...,n c , on the diagonal. Then, the closed-loop equations for joint degrees-of- 
freedom become 

w J a + [2C,fi,]t^ + = 0 

w J p + [2CA]^ + = 0 

which show that the closed-loop poles associated with joint degrees-of- freedom are stable 
poles located at ± -jfi.VT — C 2 > * = 1,2, ...,n c . The augmented joint actuator 

input commands in Eq. (27) can also be written as 


[o <5" 

f 4 1 
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l X <* J 
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(28) 
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It is necessary to show that the introduction of augmented joint actuator commands 
does not affect stability of the closed-loop system. In physical coordinates, the modified 
a-substructure state equation (i.e., x a = A am x Q + 7?' m 4) including the augmented joint 
input takes the form 


4 I 
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(29) 


in which x denotes those partitions that are in general not zero. The above form is 
obtained from identifying nonzero terms in the open-loop A a , B^, B^ y and C J a matrices 
(in physical coordinates) and by using the definitions in Eq. (20). Equation (29) clearly 
shows that the joint degrees-of- freedom are not controllable by the interior actuators, u *. 
It also shows that the augmented joint input, u given in Eq. (27) will not affect stability 
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of the closed-loop system, because of the zeros in the lower-left quarter of A am matrix. 
The closed-loop system poles are the union of eigenvalues of ( A ff m + BjJ^Gj/) and the 
closed-loop poles associated with joint degrees-of-freedom, i.e., — ± — j’0,yl — 

i = 1,2, ...,n c . However, Eq. (29) also shows that dynamics of the interior degrees-of- 
freedom are affected by the joint degrees-of-freedom, because of the nonzeros in the upper- 
right quarter of A am matrix and the nonzero in the upper part of the B J a . This interaction, 
although does not affect stability, can degrade the closed-loop system performance. In 
order to achieve satisfactory closed- loop performance, the design parameters ( and f), 
should be assigned such that vibration of joint degrees-of-freedom decays much faster 
than that of interior degrees-of-freedom. 

In summary, the procedure for designing decentralized controller gain matrix con- 
sists of three steps. First, the feedback commands for interior actuators are set to be 
u * = x q = | GjJ 0 j x a and Up = GpXp = [ Gp 1 0 ] xp with G ! J and Gp 1 designed 
based on X l = Al Jm X i + BilnU 1 * and X 0 = A 0m X 0 + B 0rr, U 0- Second, the feedback com- 
mands for the joint actuators are set to be u J a = ( G ^ + G J a ) x a and Up = (Gp + Gp) xp, 
where G J a and Gp are defined in Eq. (23) and G J a and Gp are defined in Eqs. (27) and 
(28). This set of feedback commands can be expressed by one equation: 



Finally, after the two substructures are connected together, a global control gain matrix 
for the global (or assembled) structure can be obtained by substituting Eq. (30) into 
Eq. (15). The result is 

u = Tiu = TiGx = Gx (31) 

The closed-loop poles associated with the joint degrees-of-freedom are not changed by 
the interconnection of substructures. This is proved by substituting u J = u J a + Up into 
Eq. (9), which yields y J = w J = — [2&fl, j w J — jur 7 . The set of feedback commands 
derived in this section will be called the Joint Decoupling Actuator Commands, because 
they decouple the joints from the interior degrees-of-freedom. 

IV. The Design of Decentralized Observers 

In most practical situations, the number of sensors is limited and thus the full mea- 
surement of the states is not available. In order to implement state feedback control, an 
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observer is required to reconstruct the entire system state. In this section, the procedure 
for designing decentralized observers is presented. 

The derivation is again based on the coupled substructure state-space model in 
Eq. (10), whose short notation expression is given by Eq. (11). Let the observer be 
described by 

q = Aq -\- Bu - F [y - Cq — Du) (32) 

where q = j j is the observer state vector with q a estimates x a and q 0 estimates x 0 , 

respectively. F is the observer gain matrix to be designed. According to Eq. (30), state 
feedback using estimated states is done by setting 


u = Gq = 


G ! a 

G J a + G J a 


0 

0 



(33) 


Substituting Eq. (33) into the observer equation (32) and performing algebraic manipu- 
lation, we get 



<7*1 {BiS~ l D J p J 
J + | B J „S-'Di J 
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u J 




‘ CL + DLCL 
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+ F 
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f 9a 1 
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Di J S~ 1 D J / 
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D J a J S~ l D J 0 J 



0 CL 

+ DjLGi. 


. Dys-'D J Q J \ 


rrfi 
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vt 


yi 


.y l 0< 
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(34) 


where u J = (G J a q a + G 0 q 0 ) and 


CL = Cl - Df (£>")•' Ci , £>" = D'J - D" (/)")'■ Di' 

CL = Cl - D'J d , D'L = D',‘ - (C ")' 1 Dj‘ 


The above matrices in addition to the A am , A 0m , and B 1 ^ matrices given in Eq. (20) 
completely define the modified substructure state-space models. 

Equation (34) indicates that a decentralized observer design can be achieved by 
appropriate choice of observer gain matrix. First, restrict the observer gain matrix F to 
take the form 


F = 



F J a 0 
0 Fj 



(36) 
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which is similar to the form of the controller gain matrix G in Eq. (30). Then, set the 
joint output observer gain matrices to be 


Fi =-{Bi + Fi D") (D") 
F$ = - (Bj, + Fp D 1 /) {Dj 1 ) 


Substitution of Eqs. (36) and (37) into the observer equation (34) gives 

4a = (•‘derm + ^rm^a + F 7 G 7 m + F^D^G 1 ^ 9cr — F^ J/o — Va 

qp = (-d/3m + Bp m Gp + FpCp m + FlD&Gi) qp — Fpy ! p — Fp yj 


Or, in short notation, 


4a = E a q a - [ Fl F£ ] y a 
4p — Ep qp — [ Fp Fp ] yp 


(39) 


where the definitions of E a and Ep are obvious. Therefore, by restricting the observer gain 
matrix to have a decentralized form as in Eq. (36) and by setting the joint output observer 
gain matrices as in Eq. (37), the observer equation is decoupled. The interior output 
observer gain matrices F£ and Fp still need to be designed. Define a state estimation 
error vector as 


f F a i ( 9a "1 

1 e 0 J _ l xp - qp ) 


(40) 


Then, after algebraic manipulations, subtraction of the coupled substructure state equa- 
tion (11) from the observer equation (32) yields 


^a — -dam ^a T F a C am C a 
cp = Ap m cp T Fp Cp m ep 


(41) 


The above error dynamic equation is decoupled. Any existing pole placement 
algorithm can be used to determine F£ and Fp. In order for the state reconstruction to 
be effective for feedback, F* and Fp should be designed such that eigenvalues of {A am + 
FiCam.) and ( Ap m + FpCp m ) are further to the left in the complex plane than the regulator 
poles, which, as mentioned previously, are the union of eigenvalues of (d^+S^G^ 7 ) and 
(Ap m 4- BplnGj/) and the closed- loop poles associated with joint degrees-of-freedom, i.e., 
-CiD, ±jVli\Jl -a, i = l,2,...,n c . Of course, we have made an assumption that both 
(d am ,C 7 m ) and ( Ap m ,Cp m ) are completely observable. This assumption is by no means 
a restriction. In practice, it is usually possible to arrange locations of interior sensors 
to meet this assumption. Observability of the modified substructure state-space models 
and optimization of substructure interior sensor locations are future research topics. 
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Although the derivation so far has been based on a two-component structure, 
extension of the method to multi- component structure is straightforward. The controller 
formed by the substructure observers and the Joint Decoupling Actuator Commands will 
be called the Joint Decoupling Controller. In summary, the design procedure for Joint 
Decoupling Controllers can be described by the following steps. 


(1) Determine substructure state-space models 

x» = A t x a +[Bl B >}{ U Jj) 

fyi If cl 1 
W/ i ci . 


x s + 


r d” d” 

D{ 1 D J . J 


" }l<) 
l J J l ui ] 


s = 


from finite-element modeling or from system identification of the substructures. 


(2) Derive the modified substructure state-space models 


x a = A sm x , + B[ m u[ 

v! = CL x, + dL ui 


* = 


by using Eqs. (20) and (35). 

(3) Design substructure control gains: 


(a) Design based on the modified substructure subsystem = AL x^ + 
BL u[ and then, form G[ = [ G 0 ] . 

00 Set Gj — — ( Dt ! J )~‘ (d + Di'Gl). 

Choose appropriate values of (« and fi,- for joint degrees-of-freedom and then, 
form G J a = [ 0 G JJ ], where G J t J = (D/- 7 )' 1 [ -[D?] -[2^]]. 

(c) Set u[ = G[ q„ u J s = (G J t + G J ,) q s . 

(4) Design substructure observers: 


(d) Design Fj based on (A, m ,C/ m ). 

(e) Set Ff = - (B J a + Fj D") (. D J /)~\ 

(f) Form substructure observer equation: 

q. = (A m + B' m G[ + F/C/ m + F/C"G.') q. - Fj - F? 
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(5) From global feedback gain matrix for the global structure by using the input cou- 
pling matrix 7\, i.e., using the relation u = T\u (see Eqs. (30) and (31)). 


It is emphasized that although the substructure control gain matrices need to be as- 
sembled (Step (5)) for implementation on the global (or assembled) structure, there is 
no need to assemble the substructure observers. In fact, we leave the substructure ob- 
servers unassembled so that the control implementation is decentralized. However, it 
should be pointed out that in actual implementation, and y g in the substructure 
observer equations (38) are replaced by y J , since, according to compatibility condition, 
■yJ = yJ = y J . The closed-loop system of the two-component structure, obtained by 
combining the global structure state-space model in Eq. (12), the feedback law given in 
Eq. (33), and the substructure observers in Eq. (39), and by using the input and output 
coupling matrices, is given by 


x 

q 


A 

-FT 2 C 


BTiG 

E - FT 2 DTiG 



(42) 


where E = 


E a 0 

0 Ep 

The above derivation for the substructure observers is based on the “coupled” 
state-space model in Eq. (10). Therefore, as shown in the derivation the collection of 
substructure observers together with the global feedback gain matrix (i.e., T X G ) form 
a stabilizing controller for the global structure. However, it is easy to show that each 
individual substructure observer together with substructure feedback gain matrix also 
constitute a stabilizing controller for corresponding individual substructures. To be more 
precise, define a-substructure control gain and observer gain matrices to be 



(43) 


Then, by combining the “unconnected” a-substructure model in Eq. (1), the a-observer 
in Eq. (38a) and the feedback law u r a = G ! Q q a and u J a = (G^ + G^) we get the following 
closed-loop equation 




(44) 


Using the state estimation error vector defined previously, the above equation can be 
converted to 


x a 

c a 


A am + Bi m G[ + B J a G J a -B a G a 

0 A om + FlCl m 



(45) 
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which obviously is a stable system. The closed-loop poles are the same as those of the 
a-substructure when it is connected to /3. Therefore, substructure’s closed-loop poles 
are the same before and after the connection. Besides that, each substructure’s closed- 
loop response is affected by the connection only to a slight degree. If exact states are 
available for feedback (which implies the case of full state measurement), any pair of 
connecting joint degrees-of- freedom will have the same response before and after con- 
nection, providing that before connection the initial conditions of separate joints are the 
same. Consequently, each substructure’s closed-loop responses is the same whether the 
substructures are connected or not. When estimate state is used for feedback, because 
the joint degrees-of-freedom are set to have fast decaying dynamics, the responses of 
interior degrees-of-freedom will be affected by the interconnection only slightly. This 
feature makes the proposed Joint Decoupling Controllers attractive for active control of 
space structures that are required to be connected and disconnected on a routine basis. 
No controller redesign or controller shut-off-and-turn-on is necessary before, during, and 
after the entire connecting process. 

If the dimensions of substructures are still too large for control design purposes, 
model reduction can be performed to substructure state-space models. Then, a low- 
order Joint Decoupling Controller can be designed for each substructure based on the 
reduced-order substructure model. It certainly is easier to perform model reduction at the 
substructure level than at the global structure level. Although model reduction introduces 
the so-called “spillover” problem, the closed-loop system of the global structure is stable 
as long as all the substructure closed-loop systems are stable. 

V. Examples 

A. A Mass-Spring-Damper Example 

The first example is the two-component mass-spring-damper system shown in 
Fig. 2. The two substructures a and /? are to be connected by using the rigid links 
on mass (4) and mass (5). Displacement coordinates of the global structure and the 
two substructures are designated by u>,-, tu Qi , and wpi, respectively. The outputs of the 
a-substructure are: = w Q \ and = w a ^. The outputs of the ^-substructure are: 

yp = wpz and yp = wpi. All the inputs are forces. 

The method in this paper is used to design state feedback gain matrix for the 
two substructures. Then, the substructure state feedback gain matrices are assembled to 
form a global state feedback gain matrix for the global structure. The local feedback gain 
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Figure 2: A mass-spring-damper example. 

matrix G" and G^ for the interior actuators are determined by solving the following 
optimization problems: 

to O 

Mm. */a = y o (*') t Q"x 7 + Piulful dt subj. to xl = Al^xl + B^ui 

J r oo 

o { x p) T Qp Xp 4- p(u‘)T U * dt subj. to ip = A 
The regulation cost weighting matrices were chosen to be 


ill _/ 


p™ x 0 + Ba m iip 


Qa = 


KU 0 

0 Mi 7 




K l o 

0 Mp 1 


such that the first terms in the integrals represent the sum of strain energy and kinetic 
energy corresponding to the interior degrees- of-freedom of the substructures. The control 
cost weight parameter p was set to be 0.1. The design parameters for the augmented joint 
actuator commands were chosen to be Cl 2 = 25 and 2(H = 8 such that the closed-loop 
poles corresponding to the joint degree-of-freedom are -4 ± jZ. For each substructure, 
an observer is also designed according to the proposed decentralized observer design 
procedure. The observer gain matrices are determined by using a pole placement method 
such that observer poles are 12 units to the left of regulator poles in the complex plane 
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The results are summarized in Tables 1 and 2 and Figs. 4 through 7. Table 1 lists 
the control gain and observer gain matrices of each substructure. The global control 
gain matrix for the global structure is simply an assemblage of substructure control gain 
matrices. Table 2 lists the closed-loop system poles. It is seen that each substructure has 
a pair of regulator poles -4 ± ;3 associated with the joint degree-of-freedom. The global 
structure regulator poles are the union of substructure regulator poles, except that one 
pair of — 4 ± are missing because the substructure joints have been connected to form 
the global structure. The substructure observer poles are equal to regulator poles minus 
12 as the design specification required. The global structure observer poles are the union 
of substructure observer poles. 

Assume the global structure is impacted by a unit impulse force on the joint degree- 
of-freedom at time zero. This impact force gives the joint degree-of-freedom an initial 
velocity of magnitude 2.5. Thus, the initial state of the substructure is given by x 0 = 
[0,0, 0,0, 0,0, 0,0, 0,2.5, 0,0]. In practical situations, the initial state is unknown. There- 
fore, the initial state of the observer is set to be zero, i.e., q 0 = [0> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]. 

Based on these initial conditions, the global structure open-loop and closed-loop displace- 
ment response histories at w u w 2 , and w 3 degrees-of- freedom are calculated and shown 
in Fig. 4 and Fig. 5, respectively. It is seen the proposed Joint Decoupling Controllers 
successfully suppressed the vibration of the structure. In Section IV, we mentioned that 
the Joint Decoupling Controllers also operate well at individual substructure level. This 
is demonstrated by Figs. 6 and 7. Figure 6 shows the closed-loop response histories of 
the a-substructure (for the same initial conditions) when it is disconnected from the 
^-substructure. Figure 7 compares in an expanded axis scale the closed-loop response 
of DOF1 of the global structure and the closed-loop response of DOF1 of the “discon- 
nected” a-substructure. It is evident that the system response is only slightly affected 
by the interconnection of substructures. 

B. A Three- Truss Structure Example 

The second example is the two-dimensional three-truss structure shown in Fig. 3. 
The three truss components are to be connected at the joints, at which there are collo- 
cated actuators and accelerometers. Each substructure has two interior actuators located 
at arrows designated with “a” . The initial conditions is induced by an impact force of 
magnitude 1000 applied at the lower-right corner of a in the x-direction. 

A set of Joint Decoupling Actuator Commands are designed by using the pro- 
posed procedure. The interior actuator gain matrices are determined by the LQR theory 
following the formulation in the Mass-spring-damper Example. The augmented joint ac- 
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Joint 



Figure 3: A three-truss structure example 

tuator commands are also set to be such that all joint degrees-of-freedom have the same 
closed-loop poles: — 4± j3. Observer design is not included in this example because pole 
placement calculation exceeded computer capacity. (All the computations were done on 
a Macintosch Ilci computer using MATLAB software.) Therefore, exact state feedback 
is used. 

The design results are shown in Figs. 8 through 11. Open-loop displacement re- 
sponse at the lower-right corner of or-substructure in the x-direction is shown in Fig. 8. 
The open-loop response diverges because the system has rigid-body motion. Figure 9 
shows the closed-loop displacement response histories at the lower-right corner of a- 
substructure. Figure 10 shows the closed-loop displacement histories at the lower-left 
corner of ^-substructure. Note that the responses in Fig. 10 is of order 10 -2 , which is 
1/100 of the responses in Fig. 9. Evidently, the vibration of ^-substructure is suppressed 
by its own interior actuator inputs and is not affected by the vibration of ^-substructure, 
because the interaction between a and /? is eliminated. Figure 11 shows the closed-loop 
displacement responses of the left joint degrees-of-freedom. Note that the joint responses 
have magnitudes of order 10 -3 and are suppressed in very short time. Responses of the 
7-substructure are similar and therefore are omitted. 

All the closed-loop responses were calculated by using the closed-loop equation of 
the global structure. Closed-loop responses of each “unconnected” individual substruc- 
ture for the same given initial conditions were also calculated. They turned out to be 
exactly the same as the closed-loop responses included here. As discussed at the end 
of Section IV, each substructure s response is the same whether the substructures are 
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connected or not, if exact state feedback is used. This example also demonstrates the 
computational efficiency of the proposed decentralizd design method. The global struc- 
ture has 56 DOFs (112 states). A centralized design approach will involve solving a 
Riccati equation of order 112. By dividing the global structure into three-components - 
one component has 32 DOFs, two components have 16 DOFs - the size of the control 
design problem is significantly reduced. 

VI. Concluding Remarks 

A substructure-based control design procedure using the idea of Joint Decoupling 
has been developed for decentralized control of large complex flexible structures. The 
method assumes there are collocated actuators and sensors at every joint degree-of- 
freedom. Conditions of compatibility and equilibrium at the interface between sub- 
structures are used to derive a coupled substructure state-space model for the structure 
to be controlled, which leads to modified substructure state-space models to be used 
as decentralized design basis. The design concept is summarized as follows. The con- 
trol commands for interior actuators of each substructure are chosen to be localized state 
feedback. Then, the joint actuator commands are set to cancel out all the forces acting on 
the joint degrees-of-freedom, so that interactions between substructures are eliminated. 
Finally, a set of augmented commands are added to the joint actuators to stabilize the 
joint degrees-of-freedom. The advantages of the proposed method are: (1) the design 
process is completely decentralized, (2) the controllers can function at both individual 
substructure and global structure levels, and (3) a large scale control design problem is di- 
vided into several small scale subproblems. A six-degree-of-freedom mass-spring-damper 
system and a three-truss structure are used to demonstrate the proposed method. It 
shows that the proposed controller successfully suppresses the vibration of the individual 
substructure as well as the vibration of the global structure. 
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Table 1: State feedback gain matrices obtained from decentralized control design. 


Substructure a: G a = 


. = [ n n ] 


g; 

Gi + Gt 

G T a = [ 6.6628 -9.2741 2.7544 0.1350 -2.9039 -2.1168 0 0 ] 

= [00 —3 00030] 

G J a = [ 0 0 0 0 0 0 5 1.6 ] 

Fi = -[ 0.0001 0.0200 0.2841 1.4080 0.0056 0.3446 1.9427 2.7322 ] T (10 6 ) 

F* = [ 0 0 0 0 0 0 0 -1 } T 


Substructure /?: Gp 




G ! p = [ 3.9610 -7.7589 -1.5642 -2.2097 0 0 ] 

G J P = [ 0 -6 0 -0.05 6 0.05 ] 

G J p = [ 5 1.6 0 0 0 0 ] 

F l p = -[ 1.3976 0.5154 0.0091 3.8874 5.2230 0.3449 ] T (10 4 ) 

Fp = [ 0 0 0 -1 0 0 ] T 


global structure: G = 


G'„ 0 

Gt + Gt Gi + aj 
0 GJ 
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Table 2: Poles of closed- loop systems. 


Substructure a 

Regulator Poles 

Observer Poles 

-3.0178 ±>5.6888 
-0.3371 ± j*3.8531 
— 1 .6808 ±>2.0713 
—4.0000 ±>3.0000 

-15.0178 ±>5.6888 
-12.3371 ±>3.8531 
-13.6808 ±>2.0713 
-16.0000 ±>3.0000 

Substructure /? 

Regulator Poles 

Observer Poles 

-4.4318 ±>7.1943 
— 1.2173 ± >2.7689 
—4.0000 ± >3.0000 

-16.4318 ±>7.1943 
-13.2173 ±>2.7689 
-16.0000 ±>3.0000 

global 

structure 

Regulator Poles 

Observer Poles 

—3.0178 ±>5.6888 
-0.3371 ± >3.8531 
—1.6808 ±>2.0713 
—4.4318 ± jr'7.1943 
— 1.2173 ± j'2.7689 
— 4.0000 ± >3.0000 

-15.0178 ±>5.6888 
-12.3371 ±>3.8531 
-13.6808 ±>2.0713 
-16.4318 ±>7.1943 
-13.2173 ±>2.7689 
-16.0000 ± >3.0000 
-16.0000 ±>3.0000 
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Figure 6: Closed-loop response histories of a-substructure (Mass-spring-damper exam- 
ple). 



Figure 7: Comparison of closed-loop response histories of DOF 1 (Mass-spring-damper 
example). 
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Figure 8: Open-loop response at lower-right corner of a in x-direction (Three-truss ex- 
ample). 



Figure 9: Closed-loop response at lower-right corner of a (Three-truss example). 
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Figure 10: Closed-loop response it lower-left corner of fi (Three-truss example). 



Figure 11: Closed-loop responses of left joint degrees-of-freedom (Three-truss example) 
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